Zero Landau level in folded graphene nanoribbons 



E. Prada\ P. San-Jose^, L. Brey^ 
^ Instituto de Ciencia de Materiales de Madrid (CSIC), Cantohlanco, 28049 Madrid, Spain. 
^Instituto de Estructura de la Materia (CSIC), Serrano 123, 28006 Madrid, Spain 

(Dated: September 6, 2010) 

Graphene nanoribbons can be folded into a double layer system keeping the two layers decou- 
pled. In the Quantum Hall regime folds behave as a new type of Hall bar edge. We show that 
the symmetry properties of the zero Landau level in metallic nanoribbons dictate that the zero 
energy edge states traversing a fold are perfectly transmitted onto the opposite layer. This result 
is valid irrespective of fold geometry, magnetic field strength and crystallographic orientation of 
the nanoribbon. Backscattering suppression on the N — Hall plateau is ultimately due to the 
orthogonality of forward and backward channels, much like in the Klein paradox. 



Graphene is a two-dimensional conductive mem- 
brane [Ij that can be folded into the third dimension 
without degradation to its electronic and structural prop- 
erties. Folded graphene structures are recently becom- 
ing an experimental reality [2 . Curled and folded edges 
might be much more ubiquitous than previously assumed 
[3l . Their fascinating electrostatic [H |5] and electronic 
properties [6, 7 , and the possibilities afforded by the 
third dimension, are rapidly making folded structures a 
very active graphene research subfield. Folds under a 
uniform and perpendicular magnetic field are an ideal 
system to study effectively non-homogeneous magnetic 
fields, since the flux through the nanoribbon changes sign 
across the fold. This magnetic flux inversion can happen 
over very short lengthscales, something more difficult to 
achieve in semiconductor heterostructures. 

Graphene monolayers under strong magnetic fields 
have the peculiarity of supporting a zero energy Lan- 
dau level (LL) pinned at the Dirac point, where neutral 
graphene's particle and hole bands meet. The reason 
behind the pinning has been traced back to the Index 
Theorem and the Berry phase of chiral carriers moving 
under a magnetic field . This apparently innocuous 
property has a host of profound consequences [10 . It is 
responsible for graphene's anomalous Quantum Hall ef- 
fect (QHE) predicted [TT] and observed ^ shortly after 
its discovery [9 . 

In this work we study the Quantum Hall physics in 
folded graphene nanoribbons, i.e., a double layer nanorib- 
bon obtained by folding a monolayer so that the two 
layers remain decoupled (e.g. by an insulating buffer 
[13]). We find that folds effectively become a new type 
of Hall bar edge along which the current flows in the 
Quantum Hall regime. Fold states, effectively chiral one- 
dimensional channels, are protected by topology just like 
edge states along terminated boundaries. In particular, 
we analyze the propagation of the N = {) edge state, 
or 'zero edge state' (ZES), across a fold. Such an edge 
state incoming towards the along the upper boundary of 
the top layer may either get 'transmitted' to the lower 
layer, or be 'reflected' back to the same layer, see Fig. [l] 
Both possibilities are in principle allowed by symmetry. 




FIG. 1. (Color online) Above: folded graphene nanoribbon in 
which the two layers are kept decoupled. A uniform magnetic 
field creates a magnetic flux that changes sign across the fold, 
like in a magnetic domain wall (below). The N — {) chiral 
edge state on the top layer is perfectly transmitted to the 
lower layer due to orthogonality of incoming and backscat- 
tered channels. 



Our main conclusion is that, due to orthogonality of in- 
coming and backscattering channels, the ZES is however 
always perfectly transmitted to the opposite layer, irre- 
spective of fold geometry, magnetic field, energy (within 
the zero Hall plateau) and nanoribbon crystallographic 
orientation, as long as the zero field nanoribbon is metal- 
lic. Only in the case of metallic armchair and at certain 
discrete energies close to the = 1 Landau level can 
resonant backscattering to the same layer take place. 

We consider the folded nanoribbon of width W de- 
picted in Fig. [l] with a fold of radius R. In the pla- 
nar regions of the nanoribbon, the magnetic field B is 
either parallel or antiparallel to the nanoribbon's (ori- 
ented) normal. In the folded region the magnetic field 
has an in-plane component and a component normal to 
the ribbon, B-n^ that changes sign across the fold over 
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a distance ttR. Due to graphene's two-dimensional char- 
acter, the orbital motion of the carriers in the ribbon, 
responsible for Landau levels and edge states, is only 
affected by the normal component. In the absence of 
Zeeman coupling, therefore, the fold is equivalent to a 
magnetic domain wall of width ttR. Although the spin 
degeneracy of the carriers may be lifted by the Zeeman 
coupling to the total field, spin is conserved because B 
is constant in the embedding space, and trivially factors 
out of the scattering problem. 

In semiconductor systems in the Hall regime LLs carry 
a twofold spin (s) degeneracy (neglecting Zeeman). The 
effect of boundaries on LLs is to make them dispersive. 
Graphene monolayers have an extra twofold 'valley' (r) 
LL symmetry. Atomically sharp ('terminated') bound- 
aries lift the LL valley degeneracy, forming split pairs of 
edge state branches in the dispersion relation (see Fig. 
SI in the supplementary material). These branch pairs 
run roughly parallel to each other away from the Dirac 
point, so that at any given energy each LL contributes 
with two states per edge and spin. The zero LL, however, 
is different. It splits into two divergent branches of op- 
posite energies, one particle and one hole-like, at either 
side of the Dirac point. As a consequence, it contributes 
with only half the number of edge states. 

In two decoupled graphene layers, as those resulting 
from folding a monolayer nanoribbon onto itself without 
interlayer hybridization, LLs have an extra SU(2) layer 
degree of freedom (7) that is also twofold degenerate, in 
addition to spin and valley. We can choose for example 
7 = 1 for the upper layer in Fig. [l] and 7 = —1 for 
the lower one. Although terminated edges break valley 
symmetry, they do not mix layers, which results in layer- 
degenerate edge states (along the horizontal boundaries 
m Fig. [1]). Folds, in contrast, mix the two layers but do 
not break valley symmetry due to the lack of atomically 
sharp features. Therefore, layer degeneracy is lifted in 
fold states but valley degeneracy is preserved instead. 

States in magnetic domain walls (and hence in folds) 
are carefully studied in Refs. p!4ljT6] . They are the quan- 
tum analogue of the classical snake states. In said works 
the magnetic flux profile is abrupt, which would corre- 
spond to our R = limit, and the flux interface is as- 
sumed to be infinitely long, W ^ 00. These assumptions 
allow to find analytical expressions for the energies and 
wavefunctions of the states formed around the magnetic 
domain wall. Some important features of these states are: 
(a) They are one dimensional channels confined between 
the gapped bulk of the left and right regions, with a bulk 
penetration length that depends on energy and magnetic 
length Ib = {h/\eB\y^^. (b) They are chiral, i.e. they 
propagate only in one direction, just like edge states on 
terminated edges, (c) They populate both sides of the 
domain wall (i.e. both layers). In a folded nanoribbon 
R and W are finite and solutions require a numerical ap- 
proach. However, the essential physics described above 
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FIG. 2. (Color online) Tight-binding numerical simulation 
of the charge density of the = edge state as it propa- 
gates across the nanoribbon fold (central region of each plot) 
at energies close to the Dirac point. The nanoribbon is un- 
folded flat for the purpose of these plots. The case of a two 
dimensional electron gas is shown for comparison. 



remains valid. Quasi-one-dimensional channels develop 
around the fold, but they smoothly connect to the ter- 
minated boundary states. Their penetration length is 
increased by approximately ttR. In the central region 
they propagate from top to bottom in the choice of Fig. 
[TJ and they populate both layers 7 = ±1, see vertical 
stripes on Fig. [2] These facts are independent of fold 
geometry (value of R and precise shape of fold section), 
onsite potential variations or hopping amplitude modula- 
tions, for the same topological reasons as usual Hall edge 
states along terminated edges (change of the first Chern 
number respect to vacuum). 

In the context of a transport problem in a fold struc- 
ture, incoming edge states propagate along a given layer 
towards the fold. Due to the sign change of the magnetic 
flux, the electron is then forced to change direction and 
move along the fold and onto the opposite terminated 
edges. Nothing in principle precludes that the electron 
scatters onto an arbitrary superposition of the two layers 
after traversing the fold, see lower-right panel of Fig. [2j 
We show below, however, that the transmission of the 
ZES in particular is very simple and general: interlayer 
transmission is one, since backscattering onto the same 
layer becomes impossible by orthogonality of the initial 
and final states. We first provide an analytical proof for 
zigzag and armchair boundary terminations, and then 
present a numerical study for intermediate terminations 
which confirm this result in general. 

It is well known p^THTO] that edge states along zigzag 
boundaries have a well defined (position independent) 
valley polarization, and belong to opposite valleys for 
opposite propagation directions as measured along the 
nanoribbon. Denoting the upper/lower edges by a = ±1, 
and the two layers by 7 = ±1, the valley orientation 



\to,^) of zigzag states is therefore \t±i^±i) = (1,0) and 
kTi,±i) = (O'l)- Incoming and outgoing states on the 
same layer are vahey-orthogonal. Likewise, it is straight- 
forward to solve the low energy edge states with armchair 
boundary conditions [17 . The resulting valley orienta- 
tion is Itc,^) = (za7e^^^, 1)/V2. Here e^^ = e^27r(m+i)/3^ 
where integer m is the total number of sites across the 
ribbon. This number classifies the different armchair 
nanoribbons into metallic and semiconducting regarding 
their transport in the absence of a magnetic field. If 
_ 2 (metallic armchair ribbons), the four valley orien- 
tations are pairwise orthogonal: (tq,/^/ Itq,^) = (^^^-c^/^/, 
so opposite edges on the same layer are orthogonal, like in 
the zigzag case. Therefore, since valley is preserved along 
the fold propagation, backscattering onto the same layer 
is suppressed in both cases. In semiconducting armchair 
ribbons, in contrast, there is a finite overlap between all 
edge states, and backscattering can (and does) occur. 
Extending this analysis to intermediate crystallographic 
orientations is delicate. Valley polarization is no longer 
well defined, due to the contribution of a host of surface 
states, necessary to satisfy the atomic scale details of the 
boundary condition. 

We have performed numerical simulations of the prop- 
agation of edge states on the hexagonal tight-binding 
nanoribbon of width W . The nanoribbon has parallel and 
minimal boundaries [19 , but arbitrary crystallographic 
orientation. A uniform magnetic field is applied to the 
folded nanoribbon structure by a Peierls substitution. 
We employ the recursive Green's function algorithm [20] . 
A double-sweep generalization [21 allows us to obtain 
the local charge density in real space, as seen in Fig. |2j 
Our numerical model remains valid as long as R is large 
enough so as to neglect interlayer coupling, spg orbital 
rehybridization and strains in the fold. This constraint 
is in practice not very stringent, and only requires that R 
is much greater than the interatomic distance a = 1.4 A, 
the so called continuum limit. Folded structures as the 
one studied have no strain in the continuum limit, since 
the Gaussian curvature is zero, and therefore no pseu- 
domagnetic field [U [22] arises. Regarding orbital rehy- 
bridization due to curvature, it has been shown to result 
in a coupling of mean curvature to the carrier densities, 
and the formation of electric dipoles [5, ,23^, which are 
not included in our model since they are only relevant 
for R ^ a. In any case, none of these perturbations can 
destroy fold states nor their chiral properties, which are 
protected by topology. 

In Fig. [2]we plot some low energy results for the charge 
density distribution resulting from an edge state coming 
from the top layer (from the left in the domain wall pic- 
ture) and propagating through the fold. We see that in- 
terlayer transfer of the zero energy state is perfect in the 
three ribbon terminations considered: armchair (parame- 
terized by an angle ^ = 0), zigzag {0 = 7r/6) and an inter- 
mediate crystallographic orientation {0 = 0.47r/6). The 
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FIG. 3. (Color online) Conductance versus Fermi energy 
e across the nanoribbon fold for various boundary termina- 
tions. In the N = plateau, the T = 1 observed for all 
orientations represents perfect layer-index transfer. Only at 
isolated energies in the metallic armchair case, quasi-bound 
states form and induce resonant backscattering. The N — 1 
plateau has a transmission across the fold that is dominated 
by channel interference effects, which destroys the Quantum 
Hall plateaus. 



parameters used for these simulations are: W = 81nm, 
R = 16nm, Ib = S.lnm, that corresponds to a magnetic 
field B = lOT, and energy e = lOmeV. We have nu- 
merically checked that the results are valid for arbitrary 
W,R,Ib > a. 

In Fig. [3] we present the ballistic two-terminal con- 
ductance G = Gq^^To, (where Go = 2e^/h and To, 
is the interlayer transmission of edge state a). We plot 
G as a function of the edge state energy, normalized to 
the first Landau level energy ei = y/2hvF /Ib^ for three 
different crystallographic orientations of the nanoribbon 
(^ = for armchair, 7r/6 for zigzag and O.Itt/G for an 
intermediate orientation). We use the same W and Ib 
than in the previous figure, but here R = 20.3nm. For 
B = lOT, ei = O.leV. We have included for reference 
the curve for an unfolded zigzag nanoribbon (thin dotted 
line) of the same width and with the same magnetic field. 
The ribbon with zigzag edges exhibits a perfect trans- 
mission T = 1 throughout the N = plateau (thick red 
line). Other crystallographic orientations show identi- 
cal behavior (thin green line), except for semiconducting 
armchair (not shown) that does exhibit backscattering. 
Metallic armchair ribbons (dashed blue line) may also 
exhibit backscattering at certain precise energies close 
to the N = 1 subband (transmission dips). This oc- 
curs only in setups with R greater than the magnetic 
length Ib, for which quasilocalized states can develop in 
the fold region. These are weakly hybridized with the 
ZES in the edges, and may produce sharp scattering res- 
onances at which the pseudospin and valley quickly pre- 
cess as energy crosses the quasi-bound state, see Fig. S2 
in the supplementary material. This is known as reso- 
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nant backscattering, and is a generic feature of quasi- ID 
quantum systems [24 . Resonant backscattering requires, 
however, that the induced precession significantly mod- 
ifies the ZES Dirac spinor structure in the fold, which 
only happens in the armchair case (see supplementary 
information). Backscattering on higher Hall plateaus is 
not suppressed. Since each edge supports more than one 
channel, they may mix along the fold, which leads to in- 
terchannel backscattering. Indeed, at higher plateaus, 
transmission curves are quite complex, and are domi- 
nated by channel interference effects. 

We have numerically observed that the transmission on 
the = plateau is also one in the case of small mag- 
netic fields, Ib ^ for which edge states are not well 
formed and span the whole ribbon width. The picture 
emerging from this result is that if an arbitrary clean 
nanoribbon without magnetic fiux remains metallic as 
energy approaches zero, it will have perfect ZES inter- 
layer transmission under a magnetic field when folded. 
Any such a ribbon has a level crossing of forward and 
backward channels at each Dirac point that has not been 
lifted by its boundary conditions, and is also not lifted 
by adding an arbitrary space dependent magnetic field 
to the Dirac Hamiltonian. Since forward and backward 
scattering channels remain unmixed (unavoided cross- 
ing), backscattering is zero for the ZES. This result is 
independent of the details of the fold, so in particular it 
is also valid in other geometry variations, such as non- 
orthogonal folds, steps, double folds and loops (see Figs. 
S3 and S4 in the supplementary material). 

Regarding disorder, only atomic scale defects close to 
the fold that can induce valley mixing of the ZES will 
affect the interlayer transmission. Consequently, it will 
tend to be suppressed by edge roughness. Our numeri- 
cal simulations show that near-perfect transmission, how- 
ever, persists at low energies for shallow edge roughness 
(see Fig. S5 in the supplementary material). 

In conclusion. Hall regime transport through graphene 
nanoribbon folds implements, on the N = Hall plateau, 
a coherent 7r-rotation quantum gate on the layer index. 
In other words, a ZES propagating towards the fold from 
one layer is fully transferred to the opposite layer after 
crossing the fold. This is independent of energy, magnetic 
field strength, geometric details of the fold and crystal- 
lographic orientation of the nanoribbon (as long as the 
zero field nanoribbon is metallic). Isolated departures 
from this behavior only appear at certain resonant ener- 
gies for the metallic nanoribbon with armchair edges. In 
the complementary picture of a magnetic domain wall, 
backscattering of ZESs is suppressed. This behavior is 
robust, and is a consequence of symmetry properties of 
the ZES wavefunctions, analogous to quasiparticle chiral- 
ity in the absence of magnetic field. The phenomenon is 
reminiscent of the Klein paradox in graphene monolayers, 
in which backscattering of chiral quasiparticles on valley- 
preserving defects is suppressed because the normally re- 



fiected channel is orthogonal to the incoming one. Simi- 
larly, backscattering suppression of the ZES on a domain 
wall arises because the valley spinor of backward and for- 
ward channels are orthogonal, and the Dirac Hamiltonian 
in the fold does not couple valleys. 
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Quasibound states in folded nanoribbons 



SUPPLEMENTARY MATERIAL 



Fig. |S1| : Nanoribbon bandstructure in the Hall 
regime 



Zigzag 




FIG. SI. The bands of an infinite nanoribbon with armchair 
(sohd) and zigzag (dashed) boundary orientations. The Fermi 
energy (dotted) intersects the dispersive edge state branches, 
which fixes the edge state momentum q along the boundary. 

According to Bloch's theorem, an infinite graphene 
nanoribbon generated by the repetition of a unit 
supercell of width has eigenstates that are plane 
waves of momentum q in the direction x along the 
nanoribbon. The dispersion relation in the presence 
of strong magnetic field (shown Fig. SI for forward 
propagating states) is composed of flat subbands for 
g's that correspond to states localized away from the 
boundaries (most similar to bulk Landau levels), that 
become dispersive as the states move closer to the 
boundaries. When that happens, the valley degeneracy 
of the bulk Landau levels is lifted at the same time as 
the edge state branches become dispersive. The zigzag 
and armchair boundaries do so in slightly different 
ways, although both acquire a well defined (position 
independent) valley polarization. One essential feature, 
however, is independent of the boundary details: the 
Landau level at the Dirac point (zero energy) is split 
into opposite energy branches, so that the corresponding 
edge state brances intersect the Fermi energy only once. 
This is in contrast to all the other Landau levels, which 
contribute with twice the number of edge states. The 
resulting Hall conductivity axy = Gq{2N + 1) (equal to 
the two terminal ballistic conductance G) in graphene 
monolayer nanoribbons is quantized at odd multiples 
of Go = 2e^//i, where N is the number of populated 
subbands above the Dirac point. The resulting shift at 
A/" = is the hallmark of graphene's anomalous QHE. 




FIG. S2. (Color online) Charge density around the fold for in- 
coming edge state energies slightly below the N — 1 plateau, 
where evanescent modes from the latter dominate the elec- 
tronic density. The parameters used for these simulations 
are: d)W 81nm, R 4nm, Ib 12.2nm {B 4.4T), and 
energy e = 68meV; h)W = 81nm, R = 19nm, Ib = S.lnm 
[B — lOT), and energy e = 96meV. In b) the energy lies at a 
quasi-bound state in a metallic armchair nanoribbon, causing 
resonant backscattering. In the zigzag case above, resonant 
scattering does not induce backscattering due to the valley 
degeneracy of quasi-bound states. 



Figs |S3| and |S4| : Non-orthogonal folds, and other 
geometry variations 



Experimentally, graphene nanoribbon folds may be 
produced that deviate from the fold geometry studied in 



the main text, see Fig. [S3) In particular, folds that are 
non-orthogonal to the ribbon axis are possibly the easi- 
est to produce by water jet folding techniques, and are 
better candidates to implement transport experiments, 
since it is easier to contact the two layers to different 
reservoirs. Since the precise fold geometry only affects 
the scattering wavefunction around the fold, and not the 
asymptotic edge states, the transmission is perfect, like 



in the orthogonal fold case, see Fig. S4 
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FIG. S3. (Color online) Variations of the folded nanoribbon 
geometry considered in the text that might be produced ex- 
perimentally: non-orthogonal folds (a), steps (b), double folds 
(c) and loops (d). Underneath each geometry, sketch of the 
spatial charge density of the corresponding scattering states. 
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FIG. S5. (Color online) Two-terminal conductance (normal- 
ized to the conductance quantum and averaged over 16 sam- 
ples) as a function of the edge state energy (normalized to 
the first Landau level energy) in the presence of shallow edge 
roughness. The clean nanoribbon is zigzag terminated and 
has parameters: W — 81nm, R — 20.3nm and Ib — S.lnm 
{B = lOT). Boundary disorder is parametrized in terms of 
the fraction of vacancies / along the last two atomic layers 
of the boundary (as shown schematically in the inset, which 
shows the nanoribbons unfolded flat for the purpose of these 
plots). / = 10% for the squared (red) curve, / = 25% for 
the circled (green) one and / = 50% for the triangled (blue) 
one. The case of an unfolded clean zigzag nanoribbon (black 
dashed curve) is shown for comparison. 



FIG. S4. (Color online) Computed conductance across the 
four fold variations of Fig. |S3| for system parameters like in 
Fig. 3 of the main text. On the first Landau plateau, the 
transmission is perfect, like in the orthogonal fold geometry. 
Crystallographic orientation is intermediate — O.Itt/G, but 
similar results were obtained for zigzag and armchair termi- 
nations. 



Fig. |S5| : Disorder in the boundary termination 

Transport in the presence of edge roughness is numer- 
ically computed for the folded nanoribbon of Fig. 3 with 
zigzag termination. A fraction / of individual atomic 
sites are removed from the last two boundary layers, 
modeling shallow edge roughness. In Fig. [S5]the result- 
ing average interlayer transmission is plotted for increas- 
ing fraction /. It is apparent that, especially at vanishing 
energy for which the edge state moves a certain distance 



into the bulk, the edge roughness (assumed smaller than 
Ib^ the typical edge state width) has a small effect on the 
transmission. 



